Jet Physics and Event Shape Studies at HERA 
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A review is given of the latest results on jet production and studies of event shape variables in 
deep inelastic scattering at HERA. Jet cross sections studies for inclusive jet and dijet events 
are presented and compared to next-to-leading order QCD calculations. Extraction of the 
strong coupling constant a s is discussed. 



1 Introduction 
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At HERA, studies of the hadronic final state of electron-proton deep inelastic scattering (DIS) 
ep — > eX allow to test QCD over a large range of Q 2 , Q being the virtuality of the probing 
boson. We will here focus on the latest measurements of jet production cross sections and event 
shape variables. The sensitivity of the measured observables to the value of a s will be exploited 
to determine its value. 

In the Breit frame, the production of high transverse energy events can be related to the 
gluons emitted in the hard QCD process. Most of the analyses were carried out in this frame. 
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Figure 1: Inclusive cross section measurements in the Breit frame as a function of (a) Q 2 and (b) 4>f et , 
is the jet's azimuthal angle as measured with respect to the lepton scattering plane. 
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2 Jet production in DIS 



Studies of inclusive jet production provide a powerful means of investigating QCD: the DIS 
inclusive jet cross section definition does not present the infrared-sensitivity problems related 
to the jet selection, as in the case of dijet production HI. And furthermore, the inclusive jet 



cross section is very sensitive to the value of the strong coupling constant a s , hence providing a 
powerful method of extraction of its value. 

ZEUS performed an analysis of inclusive jet production in DIS in the Breit frame i. Jets 
were found in the Breit frame with the inclusive kx clustering algorithm S, and only jet cuts in 
the Breit frame were applied. Figure |l](a) shows the inclusive jet cross section as a function of 
Q 2 . The set of selection cuts is displayed in figure ||. 

Next-to-leading order (NLO) QCD calculations are able to describe the data over the 5 
orders of magnitude within 10 — 15%. The value of a s as determined from a comparison of the 
data to the NLO QCD calculations in the high Q 2 region (Q 2 > 500 GeV 2 ) isi 



a s (M z 

For comparison, HI has found a value of@ 

a s (M z 



0.1190 ± 0.0017(stat)+° : ^ (syst) +0-0026 (th) _ 



0.1186 ± 0.0030(exp)+[| : ^( th) +0.0033 (pdf) . 



(1) 



(2) 



in similar studies of inclusive jet production. Both results are in excellent agreement with the 
world average! a s {M z ) = 0.1184 ± 0.0031. 

The same ZEUS analysis was also used to measure for the first time the azimuthal angle 
distribution of jets in the Breit frame, with respect to the lepton scattering plane (figure ||(b)). 
The shape of the distribution, of the form 
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is in very good agreement with NLO QCD, and 
elements of the hard process. 

Two possible hard energy scales are natural for 
the renormalisation scale ^jlr in NLO QCD cal- 
culations: Q and the transverse energy, Et, of 
the event. HI i investigated the effect of this 
choice in dijet production at medium Q 2 . The 
study was done in the photon-proton centre-of- 
mass frame in the kinematic region defined by 
10" 4 < x < 10~ 2 and 5<Q 2 < 100 GeV 2 . In this 
region of relatively low x and Q 2 , where higher 
order corrections are expected to be larger, the 
choice of fiR is still a matter of discussion; and 
Q is of the same order as the jets transverse en- 
ergies. 

The jets obtained with the /cy-cluster algorithm 
were selected with E^f*' > 7 GeV and E^' 2 > 
5 GeV in the pseudorapidity range — 1 < rjj^ < 2. 
Figure || shows the dijet rate as a function of 
Bjorken x in different regions of Q 2 . The NLO 
calculations with n 2 R = Q 2 are able to describe 
the data but present large scale uncertainties, 
whereas the choice [i 2 R = Q 2 + E^ results in 
smaller scale uncertainties but fails to decribe 
the distributions. 



provides a detailed test of the QCD matrix 
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Figure 2: Dijet rate as a function of Bjorken x 
in different regions of Q 2 . Data is compared with 
NLO QCD calculations with 2 different choices of 
the renormalisation scale hr. 



At higher Q 2 the particular choice of /ir has a smaller impact on the perturbative QCD 
(pQCD) predictions, which describe reasonably well the dijet datai. Also the theoretical un- 
certainties become smaller. But these jet analyses usually require a large inter-jet separation: 
either the jets have a large relative transverse energy or they are selected with a large transverse 
momentum, to avoid the effects of multi-parton emissions becoming significant. About 10% of 
the DIS sample is then classified as dijet events in these "standard" jet analyses performed at 
HERA. 

HI has investigated the minimum jet separation necessary for the NLO calculation to give 
an accurate description of dijet productions. The study was undertaken in a region of relatively 
high Q 2 : 150 < Q 2 < 35000 GeV 2 and 0.1 < y < 0J. The jets were reconstructed in the 
laboratory frame with the modified Durham algorithm Dijet events were analysed by means 

of the variable 2/2 = — w^'' 0^ * s the invariant mass of the particles clustered), where • = 
2mm{Ef, E 2 ](l — cos 0y) is the relative measure of the separation between jets i and j, and 6ij 
the angle between them. 

Figure ^ shows the measured 2/2 distribution 
normalised to the inclusive DIS cross section 
&DIS- NLO QCD predictions are in good 
agreement with the data for values of 2/2 > 
0.001, with about 1/3 of the events being 
classified as dijet events. But it overestimates 
the data dramatically in the region where 
the inter-jet separations are small. Such be- 
haviour was expected since in the small y% 
region the difference between LO and NLO 
is the largest, and also the renormalisation 
scale uncertainty and the hadronisation cor- 
rections are large, making the NLO calcula- 
tions (fixed order) unreliable. The RAPGAP 
LO Monte Carlo program is seen to be in 
good agreement with the data over all the yi 
range. 

The dijet sample defined by y2 > 0.001 was then further analysed. NLO QCD calculations were 
found to describe well the distributions of several relevant jet variables. 

3 Event shape variables 

Event shape variables are interesting observables that are sensitive to the overall topology of 
the event and therefore to higher order QCD radiation. Within the QCD framework the mean 
of an event shape F is given by 

<F>{Q) = <F >PQ C D + < F yP^-Corr.^ ^ 
0(a 2 s ) f(a 3 ,a ) 

the pQCD part being calculated at present to next-to-leading order. The second term in the 
above equation, the "power correction", refers to the non-perturbative part, necessary for the 
theory to describe the data. This non-perturbative contribution which accounts for hadronisa- 
tion corrections is expressed as a power law correction, which depends on the value of the strong 
coupling constant a s and on an empirical non-perturbative parameter cxq E9. 

The mean values as a function of Q of the thrust, broadening, jet mass and C-parameter 
where studied in the Breit frame by both the ZEUSlHI and HlH3 collaborations. A general good 




-5.5 -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 



iog 10 (y 2 ) 

Figure 3: The measured 1/2 distribution normalised 
to the inclusive DIS cross section. The data were 
compared with LO and NLO QCD calculations, and 
the RAPGAP Monte Carlo program. 



agreement was found between the ZEUS and HI results. The simultaneous 2-dimensional fit 
to all mean event shapes yielded «o ~ 0.5 ± 20% and a relatively spread range of values for 
a s (Mz), suggesting the need for inclusion of higher order terms in the theoretical calculations. 

HI has also performed fits to the differential distributions of the event shapes (see El). The 
values obtained for a s and «o were found to be inconsistent with those obtained from the fits 
to the means of the event shapes. Fits using QCD resummed calculations improve significantly 
the description of the differential distributions of the event shapes (cf. figure H"). Latest 
theoretical worklij suggests that a consistent picture can be obtained when applying the concept 
of resummed QCD calculations in a simultaneous fit to all the event shape distributions. 
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Figure 4: Differential distributions for the (a) jet mass p and (b) the thrust r t E = 1 — T: fits to the Hi data 
including QCD resummed calculations with power corrections. 
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